function diff = q_0_e(m,psi,param)

    diff = q_r(param.m_e,psi,param).*exp((1-param.alpha).*2./param.sigma^2.*int_delta_exp(m,param));

end